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Abstract Understanding the relationships between species, communities, and biodiversity are important challenges 
in conservation ecology. Current biodiversity conservation activities usually focus on species that are rare, endemic, 
distinctive, or at risk of extinction. However, empirical studies of whether such species contribute more to aspects of 
biodiversity than common species are still relatively rare. The aim of the present study was to assess the contribution 
of individual amphibian species to different facets of biodiversity, and to test whether species of conservation interest 
contribute more to taxonomic, functional, and phylogenetic diversity than do species without special conservation status. 
To answer these questions, 19 000 simulated random communities with a gradient of species richness were created by 
shuffling the regional pool of species inhabiting Emei Mountain. Differences of diversity values were then computed 
before and after removing individual species in these random communities. Our results indicated that although 
individual species contributed similarly to taxonomic diversity, their contribution to functional and phylogenetic 
diversity was more idiosyncratic. This was primarily driven by the diverse functional attributes of species and the 
differences in phylogenetic relationships among species. Additionally, species of conservation interest did not show a 
significantly higher contribution to any facet of biodiversity. Our results support the claims that the usefulness of metrics 
based only on species richness is limited. Instead, assemblages that include species with functional and phylogenetic 
diversity should be protected to maintain biodiversity. 
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1. Introduction 


Biological diversity plays an important role in supporting 
ecosystem functioning and provides key ecosystem 
services to human beings (Balvanera et al., 2006; 
Cardinale et al., 2006). Unfortunately, with the expansion 
of human populations and urbanization, human-induced 
perturbations (e.g., habitat degradation, biological 
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invasions, climate change, and pollution) have caused 
widespread loss of biodiversity at both local and global 
scales (Dunne et al., 2002; Cardinale et al., 2012). 
Given that the loss of biodiversity may induce various 
ecological and economic consequences, conservationists 
have devoted much attention to biodiversity. Most 
previous conservation studies focused primarily on 
species richness (i.e., taxonomic diversity), based on 
the assumption that species can be treated as ecological 
equivalents within communities (Perrings et al., 1995). 
In recent decades, however, a growing number of studies 
indicate that species richness represents only one measure 
of biodiversity, thus demanding more comprehensive 
approaches (Naeem et al., 2012; Pool et al., 2014). It is 
now recognized that in addition to taxonomic diversity, 
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species-specific functional traits and genetic relatedness, 
mediate the stability and resilience of ecosystems 
(Villéger et al., 2010; Rudolf et al., 2014). Therefore, 
biodiversity conservation strategies should also include 
functional diversity (1.e., the range of species traits within 
a community; Pool et al., 2014) and phylogenetic history 
(i.e., evolutionary relationships among species; Faith, 
1992) of communities. 

This is especially true for communities under 
multiple disturbances in a changing world. The rapid 
and consistent response of functional diversity can 
provide advance warning of disturbances to ecosystems. 
Moreover, phylogenetic diversity has been proposed 
as a measure to predict how evolutionary plasticity 
within communities will respond to disturbances in the 
future (Faith, 1992), and to better explain ecosystem 
productivity (Cadotte et al., 2009). Incorporating multiple 
facets of biodiversity in conservation strategies can 
thus provide unique opportunities to better understand 
community structures and dynamics (Webb et al., 2002; 
Graham and Fine, 2008), to ensure the provision of goods 
and services (Diaz et al., 2007), and to identify species of 
conservation interest (Mace, 2003; Knapp et al., 2008). 

Theoretical studies suggest that the best way to 
conserve biodiversity is to protect individual species that 
preserve different biodiversity components (Devictor 
et al., 2010). Indeed, this policy is unfeasible because 
of limited resources. Instead, current conservation 
activities have focused on species that are rare, endemic, 
distinctive, or at risk of extinction (Myers et al., 2000; 
Kier et al., 2009). However, do these species of high 
conservation value contribute more to biodiversity, and 
therefore potentially induce stronger cascading effects 
on ecosystem functioning, than other species? Or do 
individual species that contribute equally to multiple 
facets of biodiversity still relatively uncommon (but see 
Devictor et al., 2010; Pool et al., 2014), especially for 
amphibian species that play a key functional role in the 
link between aquatic and terrestrial ecosystems. 

The aim of the present study is to assess the 
contribution of individual amphibian species to different 
aspects of biodiversity. We predict that amphibian species 
will contribute in different ways to taxonomic, functional 
and phylogenetic diversity. Specifically, we first explored 
the role of individual species in taxonomic, functional 
and phylogenetic diversity, which allowed us to identify 
the ecological and evolutionary importance of each 
species. Then, we tested whether species listed by the 
IUCN as being of conservation interest (i.e. endangered, 
vulnerable and near threatened species) contribute more 
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to taxonomic, functional and phylogenetic diversity 
than species without elevated conservation status (i.e., 
species of least concern). This can help us to recognize 
the conservation value of species for different facets 
of biodiversity that are not traditionally considered in 
conservation activities. 


2. Materials and Methods 


2.1. Study area and field surveys The study area, Emei 
Mountain, is located in southwest Sichuan Province, 
China. This mountain is a series of ranges spread over 
1100km’. In total, 104 nocturnal surveys were conducted 
along 52 transects (200 m x 2 m) from April to October 
in 2013 and 2014, using the same combination of 
distance sampling and quadrat sampling methods, with 
four transects sampled per night. Distance sampling is 
considered as a positive and effective method for anurans 
(Fogarty and Vilella, 2001; Funk et al., 2003), and quadrat 
sampling is effective for stream amphibians (Dodd, 2010). 
Two persons first intensively searched the transect along 
the edge of stream by turning over logs, stones, and leaf 
litter, and hand-collected all individuals (Herbeck and 
Semlitsch, 2000; Naniwadekar and Vasudevan, 2006). Ten 
quadrats of one square meter area were then randomly 
selected within the stream, all the rocks were carefully 
removed, and individual amphibians were collected using 
a hand net (Dodd, 2010). All individuals were identified 
to species, measured (snout-vent length to the nearest 
mm), photographed, toe-clipped, and then released at the 
site of capture. 

These transects, at elevations from 476—3000m, have a 
gradient of environmental conditions. Specifically, lower 
transects are close to lentic aquatic ecosystems (e.g., 
swamps and ponds), and are dominated by evergreen 
broad-leaved forest. Transects in middle and higher 
elevations are located along slowly flowing pools and 
stream tributaries, respectively, with evergreen and 
deciduous broad-leaved forests gradually replaced by 
coniferous and mixed forests. Transects were selected at 
random, but include all suitable habitats for amphibian 
species. 


2.2. Selection of functional traits Based on published 
literatures, a set of 11 functional traits reflecting 
morphology, resource use, and life history strategy were 
selected (Trochet et al., 2014; Tsianou and Kallimanis, 
2016; Table 1). These traits described the unique 
relationship between each species and the environment, 
and importantly, reflected the ecological function of 
species in those ecosystems. Qualitative and quantitative 
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Table 1 Eleven functional traits associated with morphology, reproduction, habitat, food acquisition, and locomotion (adapted from Tsianou 


and Kallimanis, 2016). 


Functional trait Type Categories Ecological roles 
Snout vent length Continuous Centimeter Morphology 
Sexual maturity Categorical 1 year, 1-3 years, >3 years Reproduction 
Clutch position Categorical Ground, lotic habitat, lentic habitat, arboreal Reproduction 
Number of eggs per clutch Continuous Mean Reproduction 
Egg size Continuous Mean Reproduction 
Duration of breeding activity Categorical Prolonged, explosive Reproduction 
Spawn site Binary Aquatic, terrestrial Reproduction 
Life history habit Binary Aquatic, terrestrial Habitat 
Time partitioning Binary Nocturnal, diurnal Habitat 
Diet Binary Carnivorous, mainly feeding on pests, mainly feed Food acquisition 
on other insects 
Mobility mode Binary Jumper, walker, runner, climber, swimmer Locomotion 


measurements were made using digital calipers and visual 
observation; values of some traits were derived from 
Fei et al. (2009). These functional traits were snout vent 
length (SVL; cm), sexual maturity (SM; 1 year, 1-3 years, 
or >3 years), clutch position (=oviposition site) (CP; 
ground, lotic habitat, lentic habitat, or arboreal), number 
of eggs per clutch (NEC), egg size (diameter) (ES; 
mm), duration of breeding activity (DBA; prolonged, 
or explosive), spawning site (SS; aquatic, or terrestrial), 
life history habit (LHH; aquatic, or terrestrial), time 
partitioning (diel activity) (TP; nocturnal, or diurnal), diet 
[DI; carnivorous (i.e., all salamander species), mainly 
feeding on pests, mainly feeding on other insects] and 
mobility mode (mode of locomotion) (MM; jumper, 
walker, runner, climber, or swimmer). Considering the 
deficiency of data on amphibian species traits (reviewed 
in Trochet et al., 2014), selection of functional traits was 
also based on the availability of complete trait data. 


2.3. Phylogenetic analyses Sequence data were obtained 
from GenBank. Sequences were edited manually in 
BioEdit v7.0.5, aligned in Mega 6.0 using the ClustalW 
algorithm with default parameters (Tamura et al., 2013), 
and were checked by eye for ambiguous alignments. 
We used Bayesian inference (BI; performed in MrBayes 
3.1.2; Ronquist and Huelsenbeck, 2003) and Maximum- 
likelihood (ML; performed under the GTRGAMMA 
model using RAxML Web Server; Stamatakis, 2006) 
to reconstruct a mitochondrial gene tree using partial 
fragments of the 12S and 16S ribosomal genes, and the 
complete t-RNA“". The best-fitting model, GTR+G, was 
developed using ModelTest 3.7 (Posada and Crandall, 
1998). For the bayesian analysis, we ran four concurrent 
Markov Chains for five million generations, sampling 
every 100 generations . The first 25% sampled trees were 
discarded as burn in. The resulting trees were combined 


to calculate posterior probabilities for each node in a 
50% majority-rule consensus tree. Branch support in the 
maximum likelihood analysis was assessed using 1000 
nonparametric bootstrap replicates. Nodes in the trees 
were considered well supported when Bayesian posterior 
probabilities were = 0.95 and ML bootstrap support was 
= 70% (Hillis and Bull, 1993). 


2.4. Statistical analyses All sampled species were 
considered to be the regional species pool. Species 
were then classified into three categories based on their 
IUCN conservation status. EN: critically endangered or 
endangered; VU/NT: vulnerable and near threatened; and 
LC: least concern. Species in the vulnerable and near- 
threatened categories were pooled to ensure a sufficient 
sample size for statistical analysis. Since the functional 
traits data were measured on both discrete and continuous 
scales, we first used Gower’s coefficient to construct 
a dissimilarity matrix of all functional traits (Pavoine 
et al., 2009). We then used a principal coordinates 
analysis (PCoA) to build a multidimensional functional 
space following Villéger et al. (2008), to define the 
position of each species in functional space. Functional 
distances using Gower’s metric were significantly 
correlated with the Euclidean distances between species 
pairs in functional space (Mantel test, 7° = 0.776, P < 
0.001), providing evidence that functional diversity can 
effectively capture community trait diversity (Pool et 
al., 2014). Following Maire et al. (2015), the first four 
synthetic principal components of the PCoA analysis 
were retained to describe the functional space, as they 
produced a mean-squared deviation index of 0.002, 
demonstrating that these components can sufficiently 
describe interspecific distances. The four axes accounted 
for 68.83% (PC1 = 24.73%, PC2 = 20.53%, PC3 = 
12.47% and PC4 = 11.10%; Figure 1 and Table 2) of the 
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Figure 1 Distribution of all species in the functional space defined by the regional species pool, based on the four first PCA axes. (A) PC1 
and PC2 of the functional space, (B) PC3 and PC4 of the functional space. EN (critically endangered or endangered), VU/NT (vulnerable 
and near threatened) and LC (least concerned) species were plotted in red, blue, and green, respectively. Colored surfaces correspond to the 
functional space occupied by each category, and crosses to the center of gravity of the convex hulls. Species codes are in Table 3. 
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Table 2 Pearson correlation coefficients between the first four synthetic principal components of the PCoA analysis axes and the 11 
functional traits. Statistically significant correlation coefficients are in bold. 


PCoA1 PCoA2 PCoA3 PCoA4 
Snout vent length 0.62 0.42 0.14 0.19 
Sexual maturity -0.26 -0.25 -0.3 -0.64 
Clutch position 0.15 -0.56 -0.25 0.43 
Number of eggs per clutch -0.09 0.87 -0.1 0.18 
Egg size 0.66 -0.49 -0.24 -0.19 
Duration of breeding activity 0.01 0.71 0 -0.21 
Spawn site -0.44 -0.24 0.04 0.65 
Life history habit -0.8 0.18 -0.04 0.03 
Time partitioning 0.12 -0.25 0.82 -0.01 
Diet -0.61 -0.21 0.58 -0.23 
Mobility mode 0.76 0.12 0.34 -0.06 


initial inertia in trait values. Each species had a distinctive 
position in functional space, which was used below to 
calculate the contribution of individual species to overall 
aspects of biodiversity. 


Contribution of individual species to multiple facets 
of biodiversity Taxonomic diversity was represented 
by the number of species, functional diversity by 
functional richness index, and phylogenetic diversity by 
the total branch length of species within the community 
(Pool et al., 2014). We simulated randomly assembled 
communities to identify the contribution of individual 
amphibian species to different facets of biodiversity 
because (1) it was easy to obtain a large number of 
communities, which produced a gradient in species 


richness; (2) the local stochastic processes involved 
in community assembly can be mimicked by random 
sampling from a regional species pool (Loreau, 2001); 
and (3) simulations allow for control of concurrent 
changes in species richness and random extinctions 
of species (Larsen et al., 2005). Based on the four- 
dimension functional space, random communities should 
contain more than four species to allow calculation of 
functional richness. Therefore, five species were first 
selected randomly from the regional species pool using 
a bootstrap procedure with 1000 repetitions (i.e., a set of 
1000 random communities, each with five species, was 
obtained). Since we assumed that the regional species 
pool contained S species in total, the same protocols 
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were then performed to obtain a set of 1000 random 
communities with six species, a set of 1000 random 
communities with seven species, and so on, up to a set of 
1000 random communities with (S — 1) species. In total, 
(S—5) x 1000 random communities were simulated. The 
abundance of each species was not considered as it did 
not affect the taxonomic, functional, and phylogenetic 
diversity indices we used. 

The contribution of each species to different aspects 
of biodiversity was quantified as follows: the taxonomic, 
functional, and phylogenetic diversity of each random 
community were calculated firstly, and then all the 
diversity indices were recalculated by removing each 
species. The differences in diversity values before and 
after species removal were computed. The relative 
contribution of each species was assessed by calculating 
the mean percentage of diversity change within each 
random communities. 


Contribution of conservation status to facets 
biodiversity The mean relative contribution to 
biodiversity of species of the same conservation status 
was calculated to explore whether species of conservation 
interest (i.e., EN and VU/NT species) contribute more 
to taxonomic, functional and phylogenetic diversity than 
species without elevated conservation status (i.e., LC 
species) do. Since the three groups have different number 


of species, we used Nemenyi—Damico—Wolfe—Dunn tests 
(Nemenyi, 1963) (a nonparametric Tukey-type, multiple 
comparison post hoc test) to compare the difference of 
biodiversity contribution between pairwise categories, 
All statistical analyses were conducted in R 3.3.2 (R 
development Core Team, 2017). 


3. Results 


A total of 24 amphibian species were detected in the 
regional species pool, with 5 EN species, 5 VU/NT 
species, and 14 LC species. The number of individual 
animals of each species ranged from 1 to 10, with a mean 
of 4.2 + 3.2 SD. The most abundant amphibian species 


were Spiny-bellied frog (Quasipaa boulengeri) and 
Baoxing treefrog (Rhacophorus dugritei), which made 
up 31.75% of the total individuals. The rarest species 
was Longdong stream salamander (Batrachuperus 
londongensis), with only one individual found. 


3.1. Contribution of individual species to biodiversity 
All 24 species contributed similarly to taxonomic 
diversity, but differently to functional and phylogenetic 
diversity of simulated random communities (Table 3). 
Specifically, the Asiatic toad (Bufo gargarizans), the 
Omei music frog (Nidirana daunchina) and the Spiny- 
bellied frog (Quasipaa boulengeri) were the species 


Table 3 Contribution of individual species to different components of biodiversity (TD: taxonomic diversity; FD: functional diversity; PD: 


phylogenetic diversity). 


Scientific name Code IUCN level TD (%) + SE FD (%) + SE PD (%) + SE 
Batrachuperus londongensis bal EN 6.70 + 0.02 11.66 + 0.08 7.06 + 0.04 
Oreolalax omeimontis oro EN 6.44 + 0.02 5.76 + 0.07 10.34 + 0.05 
Quasipaa boulengeri qub EN 6.77 + 0.02 16.85 + 0.09 7.79 + 0.04 
Scutiger chintingensis scc EN 6.18 + 0.01 1.57 + 0.04 6.06 + 0.04 
Vibrissaphora boringii vib EN 6.57 + 0.02 1.47 + 0.04 6.88 + 0.04 
Batrachuperus pinchonii bap VU/NT 5.94 + 0.01 7.18 + 0.08 7.69 + 0.04 
Megophrys omeimontis meo VU/NT 6.59 + 0.02 8.15 + 0.06 6.86 + 0.04 
Oreolalax major orma VU/NT 5.94+ 0.01 2.15 + 0.03 7.39 + 0.04 
Oreolalax multipunctatus ormu VU/NT 6.82 + 0.02 11.28 + 0.09 9.45 + 0.05 
Oreolalax schmidti ors VU/NT 6.72 + 0.02 6.11 + 0.07 8.72 + 0.05 
Amolops granulosus amg LC 6.37 + 0.02 3.78 + 0.06 10.83 + 0.06 
Bufo gargarizans bug LC 6.95 + 0.02 35.76 + 0.06 7.17 + 0.04 
Hyla chuanxiensis hyc LC 6.23 + 0.01 9.47 + 0.07 9.52 + 0.04 
Megophrys shapingensis mes LC 6.32 + 0.02 1.09 + 0.03 7.45 + 0.05 
Nidirana adenopleura nia LC 6.04 + 0.02 9.08 + 0.08 9.62 + 0.06 
Nidirana daunchina nid LC 6.47 + 0.02 21.27+0.11 9.78 + 0.06 
Odorrana schmackeri ods LC 6.86 + 0.02 4.28 + 0.06 9.57 + 0.05 
Odorrana margaretae odm LC 6.26+0.01 11.41 + 0.09 10.26 + 0.06 
Oreolalax popei orp LC 6.42 + 0.02 9.56 + 0.07 9.18 + 0.05 
Paramegophrys oshanensis pao LC 6.24 + 0.02 12.49 + 0.08 6.47 + 0.04 
Rana omeimontis rao LC 6.38 + 0.02 3.14 0.06 10.23 + 0.06 
Rhacophorus chenfui rhe LC 6.78 + 0.02 6.03 + 0.06 8.42 + 0.04 
Rhacophorus dugritei thd LC 6.47 + 0.01 11.10 + 0.08 8.12 + 0.05 
Rhacophorus omeimontis tho LC 6.65 + 0.02 4.51 + 0.06 8.15 + 0.04 
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Figure 2 Box plots of each IUCN conservation status’ (EN: 
critically endangered or endangered; VU/NT: vulnerable and 
near threatened; and LC: least concern) contribution to different 
components of biodiversity (TD: taxonomic diversity; FD: 
functional diversity; and PD: phylogenetic diversity). Error bars 
indicate the standard error. 


that contributed the most to functional diversity (i.e. 
35.76% + 0.16% SE, 21.27% + 0.11% SE, and 16.85% 
+ 0.09% SE; respectively). Species that contributed the 
least to functional diversity were the Shaping horned toad 
(Megophrys shapingensis; 1.09% + 0.03% SE), the Omei 
mustache toad (Vibrissaphora boringii; 1.47% + 0.04% 
SE) and the Jinding odontoid toad (Scutiger chintingensis; 
1.57% + 0.04% SE). In terms of phylogenetic diversity, 
the Sichuan sucker frog (Amolops granulosus), the 
Omei tooth toad (Oreolalax omeimontis) and the Green 


odorous frog (Odorrana margaretae) were the substantial 
contributors (i.e. 10.83% + 0.06% SE, 10.34% + 0.05% 
SE, 10.26% + 0.06% SE; respectively), while the Jinding 
odontoid toad (Scutiger chintingensis), the Oshan 


metacarpal-tubercled toad (Paramegophrys oshanensis) 
and the Omei horned toad (Megophrys omeimontis) 
contributed less (6.06% + 0.04% SE, 6.47% + 0.04% 
SE, 6.86% + 0.04% SE; respectively) than other species 
(Table 3). 


3.2. Contribution of conservation status to biodiversity 
The relative contributions of EN species to taxonomic, 
functional and phylogenetic diversity were 6.52% + 
0.04% SE, 7.46% + 1.20% SE, 7.63% + 0.29% SE; 
respectively. The relative contributions of VU/NT 
species to taxonomic, functional and phylogenetic 
diversity were 6.40% + 0.08% SE, 6.69% + 0.59% 
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SE, 8.02% + 0.19% SE; respectively. For LC species, 
the relative contributions to taxonomic, functional and 
phylogenetic diversity were 6.43% + 0.02% SE, 10.21% 
+ 0.62% SE, 8.70% + 0.10% SE; respectively. Nemenyi-— 
Damico—Wolfe—Dunn multiple comparison post hoc 


tests indicated that there were no significant differences 
between pairwise categories of conservation status of the 
relative contribution to taxonomic diversity (EN - VU/ 
NT: P = 0.88, EN - LC: P = 0.76, VU/NT - LC: P = 0.94; 
respectively), functional diversity (EN - VU/NT: P = 1.00, 
EN - LC: P = 0.83, VU/NT - LC: P = 0.86; respectively) 
and phylogenetic diversity (EN - VU/NT: P = 0.94, EN 
- LC: P = 0.31, VU/NT - LC: P = 0.48; respectively) 
(Figure 2). 


4. Discussion 


Our results revealed that amphibian species contributed 
similarly to taxonomic diversity, but differently 
to functional and phylogenetic diversity. This was 
primarily driven by the diverse functional attributes of 
individual species present in the functional space, and the 
phylogenetic relationships between species. 

Traditional conservation theories suggested that the 
maximum protection of taxonomic diversity (species 
richness) can preserve multiple facets of biodiversity 
simultaneously (Devictor et al., 2010). However, it is 
impossible to protect all the species simultaneously, and 
it has been suggested that endemic and distinctive species 
be given conservation priority (Myers et al., 2000; Kier 
et al., 2009). Given that not all endemic and distinctive 
species contributed substantially to taxonomic, functional 
and phylogenetic diversity in the present study, we argue 
that effective conservation strategies should not only 
focus on species per se, but that such strategies need a 
comprehensive method that includes species richness, 
functional traits, and their evolutionary relationships 
(Brooks, 2006; Lee and Jetz, 2008). 

In fact, one of the important research areas in 
conservation ecology was the unraveling of linkages 
between species traits, genetic relatedness, and extinction 
risk, which were critical criteria for prioritizing 
conservation efforts (Murray and Hose, 2005). Our 
results showed that amphibian species that contributed 
the most to functional diversity usually exhibited specific 
functional attributes (e.g., the Omei music frog and 
Spiny-bellied frog), demonstrating that these species may 
play distinct functional roles in ecosystems. These species 
displayed larger body size, smaller number of eggs per 
clutch, larger egg size, and shorter breeding period. More 
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importantly, these traits may lead to higher vulnerability 
when species are under various types of disturbances. 
This can be explained by the fact that species with larger 
body size generally have a longer life cycle (Cardillo, 
2005), and species with low fecundity have a slower 
population turnover, thus impeding rapid population 
restoration after historical contractions or catastrophic 
population crashes (Williams and Hero, 1998; Hero et 
al., 2005; Sodhi et al., 2008). Moreover, these species are 
mainly carnivorous, with poor capacity for locomotion, 
leading to high predation risk and over-exploitation by 
humans (Collins and Crump, 2009; Chan et al., 2014). In 
addition, our findings were also consistent with Lips et al. 
(2003), showing that riparian amphibians were generally 
more easily affected than terrestrial species, probably 
because many potential causal factors of population 
decline are stream-borne, such as disease and pollution 
(Williams and Hero, 1998). Similarly, species that 
contributed the most to phylogenetic diversity usually 
showed distinctive evolutionary relationships (e.g., 
the Sichuan sucker frog and Omei tooth toad). These 
“phylogenetically rare” species should be considered 
to have priority for protection as their extinction will 
cause larger loss of evolutionary information within the 
community (Winter et al., 2013). Interestingly, we found 
that the degree of relative contribution of a species to 
functional and phylogenetic diversity were not similar. 
This result indicated that species within a community 
may have redundant functional traits, but show a unique 
position in the phylogeny (e.g. Amolops granulosus and 
Rana omeimontis), or species with distinct functional 
traits may represent phylogenetic conservatism within a 
community (e.g., Bufo gargarizans and Paramegophrys 
oshanensis; Pool et al., 2014). All of these observations 
demonstrated that the conservation importance of species 
was affected by both their functional and phylogenetic 
attributes (Faith, 1992; Diaz and Cabido, 2001). And 
importantly, our results may give potential conservation 
priority for each species, with higher priority given to the 
conservation of species contributing most to functional 
and/or phylogenetic diversity in a given region. 

Species in categories of conservation interest did 
not display a significantly higher contribution to any 
aspect of biodiversity in our study area. Indeed, the local 
extirpation of either endangered or least concern species 
can result in dramatic loss in specific facet of community 
biodiversity. For instance, Quasipaa boulengeri (EN) 
and Bufo gargarizans (LC) contributed substantially to 
functional diversity, while Oreolalax omeimontis (EN) 
and Amolops granulosus (LC) contributed substantially 


to phylogenetic diversity. Therefore, a set of species 
assemblages that include diverse threatened and common 
species should be protected to maintain a variety of 
biodiversity facets (Pool et al., 2014). 

In conclusion, the present study revealed the 
contribution of individual species to multiple biodiversity 
facets. Our results suggested that species of conservation 
interest were not always the greatest contributors to 
taxonomic, functional, and phylogenetic diversity. 
Moreover, our results suggest that the integrated diversity 
components of threatened and common species within a 
community may help preserve community biodiversity. 
In turn this can allow conservationists to optimize 
conservation strategies and protect biodiversity effectively 
(de Carvalho and Tejerina-Garro, 2015). Given that this 
study was constrained to a limited area (1.e., local scale), 
and biodiversity crisis is a global problem in space/time, 
how to maintain biodiversity across larger temporal 
and spatial scales in a changing world should be further 
studied (Gering et al., 2003). 
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